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The  stability  of  Pt-skin  on  Pt3Co  alloy  against  Co  and  Pt  dissolution  with  respect  to  Pt-skin  thickness, 
surface  structure  and  composition  is  investigated  using  density  functional  theory  calculations.  It  is  found 
that,  even  under  oxygen  environment,  Co  migration  to  the  surface  is  suppressed  by  the  thicker  Pt-skin  (2 
-3  Pt  atomic  layers),  confirming  experiments.  However,  the  instability  of  single  Pt  layer  is  attributed  to 
the  electronic  effect  of  oxygen.  The  adsorbed  oxygen  redistributes  charge  from  Pt— Co  region  to  Pt— 0 
region,  weakening  the  Pt— Co  bond  and  facilitating  Co  migration  to  surface  via  Co— 0  bond  formation.  We 
further  note  that  in  this  system,  the  Co  migrates  easier  in  (111)  surface  than  in  (100),  attributed  to  the 
difference  in  the  oxygen  adsorption  sites  (fee  vs  bridge).  A  minimal  negative  electrode  potential  shift  of 
0.06-0.09  V  for  Pt  dissolution  is  noted  for  thicker  Pt-skin  systems,  indicating  stability  close  to  pure  Pt. 
The  Pt-skin  composition  is  varied  by  introducing  different  Migration  Barrier  Layers  (MBLs)  =  Ru,  Rh,  Pd, 
Os,  Ir,  in  the  mid-Pt-skin  region  and  found  that  when  MBL  is  Os  or  Pd,  a  novel  low-cost  composition  and 
more  stable  MBL-substituted  Pt-skin/Pt3Co  ORR  catalyst  emerges. 

©  2013  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Overall,  the  cost,  performance  and  durability  of  the  catalyst  are 
the  three  most  important  factors  that  affect  commercialization  of 
polymer  membrane  electrolyte  fuel  cell  (PEMFC).  It  is  well-known 


*  Corresponding  author. 

E-mail  address:  mcescano@u-fukui.ac.jp  (M.C.S.  Escano). 


that  the  kinetics  of  the  oxygen-reduction-reaction  (ORR)  at  the 
cathode  is  still  sluggish  even  on  Pt.  At  present,  the  most  promising 
near-term  electrocatalyst  for  the  ORR  are  seen  to  be  Pt-based  alloy 
systems  [1-4].  Among  these,  the  Pt-Co  systems  have  received  a 
considerable  attention  because  of  high  ORR  activity  [2,4]  and  better 
stability  in  acidic  environment  [5-7].  Thus,  further  nanostructuring 
of  Pt-Co  systems  for  an  enhanced  performance  and  stability  as 
compared  to  the  conventional  pure  Pt  catalysts,  while  maintaining 
a  low  Pt  loading,  is  of  tremendous  interest.  However,  the  loss  of  the 
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3d  transition  metals  through  easy  formation  of  surface  oxides  is 
known  to  significantly  reduce  the  ORR  activity  in  Pt  alloys  after 
some  potential  cycling,  and  therefore  limits  the  use  of  these  sys¬ 
tems  in  PEMFC  [1-4].  It  has  been  observed  experimentally  that  a 
typical  pure  Pt  outerlayer,  often  called  “Pt-skin”,  is  formed  on  Pt— 
Co  systems  after  electrochemical  stabilization  of  the  As-prepared 
alloy  or  under  electrochemical  measurements  [1,4, 8, 9].  Such  Pt 
layer  formation  on  the  alloy  surface  was  also  found  in  other  Pt-3d 
transition  metal  combinations  [3,4,10,11  .  Despite  the  fact  that  the 
Pt-skin  (“outerlayer”  and  “skin”  are  used  interchangeably)  can 
serve  as  a  protective  covering,  Co  dissolution  is  not  significantly 
eliminated  [8  .  It  was  shown  by  STM  image  that  the  Pt  outerlayer 
has  (111)  surface  structure  8].  However,  more  recent  experiment, 
involving  the  same  (111 )  surface,  but  with  a  thicker  Pt-skin  (i.e.  2  Pt 
layers)  formed  by  nanocapsule  method,  reports  better  stability 
against  Co  dissolution  12].  Next,  another  experiment  reports  su¬ 
perior  stability  against  Co  dissolution  upon  potential  cycling  for 
annealed  nanoparticles  of  5-7  nm  sizes.  These  systems  exhibit  a 
(100)  surface  structure  of  the  Pt-skin  and  is  composed  of  2-3 
atomic  Pt  layers  [13].  Moreover,  the  loss  of  electrochemical  active 
surface  area  (ECSA)  is  found  to  be  minimal  even  after  5000  cycles. 
These  experimental  researches  have  motivated  the  present  study  to 
explore  the  atomic  level  properties  and  energetics  of  Co  dissolution 
with  respect  to  the:  (1 )  Pt-skin  thickness  and  (2)  surface  structure, 
which  has  not  been  explored.  A  first-principles  study  based  on 
density  functional  theory  (DFT)  is  conducted  on  the  migration 
behavior  of  Co  from  the  alloy  substrate  to  the  surface  of  the  Pt-skin, 
accompanied  by  Co-0  bond  formation  on  surface.  Based  on  the 
experiments  on  Pt-Co  systems  [8,12,13],  alloy  surfaces  are  often 
found  to  be  decorated  by  Co  oxides,  which  eventually  dissolve  into 
the  acid  solution.  Theoretical  works  have  explored  the  dissolution 
of  3d  transition  metals,  mostly  Co,  Ni  and  Fe,  in  Pt-3d  alloys  as  well 
as  in  the  Pt  monolayer  on  Pt-3d  alloys  [14-16].  These  studies  have 
reported  ease  of  3d  transition  metal  migration  with  O  coverage  as 
compared  to  that  of  the  zero  coverage.  More  recently,  the  oxygen 
adsorbate-induced  segregation  of  Co  has  been  reported  for  Pt-Co 
alloy  at  nanoscale  [17].  Since  the  mechanism  for  this  observation  is 
also  lacking,  it  is  therefore  necessary  to  seek  the  origin  of  the  in¬ 
fluence  of  oxygen  in  the  dissolution,  in  particular,  the  Co,  within  an 
electronic  perspective.  For  completeness  of  stability  evaluation,  the 
electrochemical  potential  shift  of  Pt  dissolution  in  the  thicker  Pt 
skin  systems  and  of  different  surface  structure  is  further  evaluated 
and  compared  with  the  pure  Pt  catalysts.  Lastly,  a  quest  for  new 
design  of  Pt-skin  has  prompted  this  research  to  undertake  a  non-Pt 
mid-layer  substitution  in  the  thicker  Pt-skin  systems.  The  mid-layer 
(herein  termed  as  Migration  Barrier  Layer  or  MBL)  is  composed  of 
either  one  of  the  following  metals:  Ru,  Rh,  Pd,  Os,  Ir.  Basically,  these 
4d  and  5d-transition  metals  have  been  considered  in  magnetic 
multilayers  in  spintronics  as  “capping”  layers  to  prevent  oxidation 
or  intermixing  of  metals  [18-20].  Here,  the  reactivity  and  stability 
against  Co  and  Pt  dissolution  of  the  MBL-substituted  Pt-skin  cata¬ 
lysts  under  oxygen  environment  are  evaluated.  This  screening  of 
MBLs  could  give  way  for  a  viable  Pt-skin  composition  on  Pt3Co 
alloy.  It  is  viewed  that  this  present  study  provides  new  insights  into 
the  migration  of  Co  to  the  surface  with  respect  to  different  Pt-skin 
structures,  thickness  and  composition,  thereby  aid  in  furthering  the 
design  of  stable  Pt-skin/Pt3Co  ORR  catalyst  for  PEMFC  cathode.  Also, 
a  novel  and  potentially  more  stable  and  more  reactive  Pt-skin 
composition  on  Pt-Co  alloy  is  drawn. 

2.  Computational  methods 

The  first-principles  study  is  conducted  within  the  spin- 
polarized  density  functional  theory  (DFT)  [21,22]  framework  us¬ 
ing  Projector  Augmented  Wave  (PAW)  method  [23]  and  plane-wave 


basis  set.  The  calculations  are  performed  using  the  spin-polarized 
version  of  Vienna  ab  initio  Simulation  Package  (VASP)  [24-27] 
which  can  properly  treat  magnetic  character  of  O2  and  Co.  The 
Pt-Co  alloy  surface  is  modeled  using  Ll2-ordered  Pt3Co  bulk 
structure  cut  in  (111)  and  (100)  directions.  The  Ll2  structure  of  the 
bulk  alloy  is  composed  of  Pt  atoms  at  the  face  centers  and  the  Co 
atoms  at  the  corners  of  the  cell.  Therefore  two  slabs  have  been 
formed:  Pt3Co(lll)  and  Pt3Co(100).  The  Pt-skin  on  Pt-Co  systems 
are  modeled  by  laying  1,  2,  or  3  atomic  layers  of  Pt  pseudo- 
morphically  on  Pt3Co(lll)  and  Pt3Co(100)  slabs.  This  enables 
attainment  of  the  same  surface  structure  as  the  slab.  The  multiple 
layers  of  the  Pt-skin  also  attain  a  fee  stacking.  The  surface  unit  cell 
size  is  (2  x  2)  and  the  vacuum  size  is  -12.0  A.  The  optimized 
surfaces  are  obtained  by  relaxing  the  Pt  layers  and  the  two  top¬ 
most  layers  of  the  Pt3Co  substrates.  The  optimized  structures  are 
shown  in  Figs.  1  and  2  (for  side  views  only,  the  corresponding  top 
views  are  given  in  Figs.  5  and  7)  and  are  denoted  as  either  PtnML/ 
Pt3Co(lll)  and  PtnML/Pt3Co(100)  (depending  on  the  facet  of  the 
Pt3Co  substrate),  where  n  is  the  number  of  Pt  atomic  layers 
composing  the  Pt-skin.  Details  of  the  atomic  structure  and  stability 
evaluation  are  discussed  in  the  next  section.  At  this  point  however, 
it  worthwhile  to  mention  that  the  model  systems  {PtnML/Pt3Co(lll ) 
and  PtnML/Pt3Co(100)}  appropriately  describe  the  multiple  Pt-skin 
layers,  the  resulting  surface  structures  and  the  uniform  Pt3Co 
substrate/core  composition  observed  in  experiments  [12,13  as 
previously  discussed  in  the  introduction.  Next,  A  plane-wave 
expansion  with  a  cutoff  energy  of  400  eV  is  used  throughout  the 
calculations.  The  conjugate-gradient  method  [28]  is  used  for  the 
relaxations  and  the  positions  of  the  atoms  are  moved  until  the 
maximum  force  (on  each  atom)  is  converged  to  less  than 
0.02  eV  A-1.  The  generalized  gradient  approximation  (GGA)  based 
on  the  Perdew,  Burke  and  Ernzerhof  (PBE)  functional  [29]  is  used 
for  the  exchange-correlation  correction.  Brillouin-zone  integrations 
are  performed  on  a  grid  of  (4  x  4  x  1 )  Monkhorst-Pack  k  points 
with  a  smearing  of  Methfessel— Paxton  method  [30].  These  calcu¬ 
lation  parameters  are  tested  and  rendered  an  excellent  agreement 
of  equilibrium  lattice  parameters  of  fcc-Pt,  fcc-Co  and  Ll2-Pt3Co  {Pt: 
3.97  A  Co:  3.52  A,  Pt3Co:  3.89  A}  with  experiments  {Pt:  3.96  A  31  ], 
Co  3.55  A  [32],  Pt3Co:3.85  A  [33]}.  For  the  oxygen,  the  DFT-GGA 
implementation  yields  a  bond  length  of  1.24  A  and  a  dissociation 
energy  of  5.67  eV,  in  general  agreement  with  experiment  (1.21  A; 
5.25  eV)  [34]  and  with  other  DFT-GGA  calculations  (1.24  A; 
5.64  eV)  [35,36]. 

3.  Results  and  discussions 

3.1.  Atomic  structures  of  PtnMjPt3Co  systems 

As  discussed  previously,  the  Pt-skin  on  Pt-Co  alloys  are  repre¬ 
sented  by  PtnML/Pt3Co  systems,  where  n  =  1, 2,  or  3  Pt  atomic  layers. 
Thus  a  total  of  six  PtnML/Pt3Co  structures  are  considered  in  this 
study  -  three  structures  belonging  to  each  of  (111 )  and  (100)  facets. 
The  systems  are  denoted  as:  PtiML/Pt3Co(lll),  Pt2ML/Pt3Co(lll), 
and  Pt3ML/Pt3Co(lll)  for  the  (111)  facet  as  shown  in  Fig.  1(a)— (c) 
and  PtiML/Pt3Co(100),  Pt2ML/Pt3Co(100),  and  Pt3ML/Pt3Co(100)  for 
(100)  facet  as  shown  in  Fig.  2(a)— (c).  The  pure  Pt  outerlayers  are 
labeled  as  1,  2  and  3  beginning  from  the  top  of  the  Pt3Co  substrate. 
Regarding  the  Pt3Co(lll)  substrate,  each  layer  is  composed  of  3  Pt 
atoms  and  1  Co  atom  as  shown  in  Fig.  1.  The  relevant  bond  lengths 
and  interlayer  distances  are  given  from  the  2nd  to  the  5th  column 
of  Table  1.  The  first  interlayer  distance  (do-i)  are  generally  larger 
than  that  of  pure  Pt(lll).  The  rest  of  the  interlayer  distances  (di_2i 
d2_ 3)  show  a  quite  wider  range  of  values  across  different  Pt-skin 
thickness.  However,  Pt-Co  bond  length  (dPt_co)»  remains  un¬ 
changed  with  respect  to  the  Pt-skin  thickness.  The  lateral  Pt-Pt 
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Pt3ML/Pt3Co(lll) 


Fig.l.  Optimized  atomic  structure  of  PtnML/Pt3Co(lll).  (a),  (b)  and  (c)  correspond  to  n  =  1,  n  =  2  and  n  =  3,  respectively.  The  Pt3Co(lll)  substrate  is  indicated  in  (a).  The  interlayer 
distances,  d0-i,  di_2  and  d2_ 3,  the  Pt  surface,  Pt  mid-layers  and  the  dPt-co  bond  are  shown. 


distance  in  the  alloy  system  (2.75  A)  is  shorter  than  that  of  the  pure 
Pt  (2.81  A).  In  the  case  of  the  PtnML/Pt3Co(100)  structures  {Fig.  2(a)— 
(c)},  almost  the  same  trend  in  the  interlayer  distances  can  be 
observed  (7th  to  the  11th  column  of  Table  1).  The  do-i  is  also  larger 
than  that  of  pure  Pt(100)  and  there  is  also  wide  range  of  values  for 


the  upper  interlayers  (di_2,  CI2-3).  However,  contrary  to  PtnML/ 
Pt3Co(lll),  the  dpt-co  varies  with  the  thickness  of  the  Pt-skin.  The 
dpt-co  is  larger  for  thicker  Pt-skin.  The  lateral  Pt-Pt  distance  is  the 
same  as  that  in  pure  Pt(100),  2.81  A  It  should  be  noted  that  the 
Pt3Co(100)  substrate  is  obtained  by  interchanging  the  Co  and  Pt 


(lOO)-cut  Ll2  Pt3Co 

Pt3ML/Pt3Co(100) 


Pt2ML/Pt3Co(100) 


Pt1ML/Pt3Co(100) 


Pt3Co(100) 

substrate 


surface 


(a) 


surface 


Fig.  2.  Optimized  geometry  of  PtnML/Pt3Co(100).  (a),  (b)  and  (c)  correspond  to  n  =  1,  n  =  2  and  n  =  3,  respectively.  The  Pt3Co(100)  substrate  is  labeled  in  (a).  The  interlayer  distances, 
do-i,  di_2  and  d2_3,  the  Pt  surface,  Pt  mid-layers  and  the  dPt-co  bond  are  shown.  The  inset  figure  shows  the  (lOO)-cut  Ll2-Pt3Co  with  Co  shift  in  position  to  attain  consistency  of  Co 
composition  per  layer  as  that  used  in  Pt3Co(lll)  substrate. 
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Table  1 

Atomic  structure  {i.e.  interlayer  distances  (d0-i,  di_2,  d2- 3),  Pt-Co  bond  length  (dPt-co)  and  lateral  Pt-Pt  distances}  of  the  PtnML/Pt3Co(lll),  PtnML/Pt3Co(100),  Pt(lll)  and 
Pt(100).  Experimental  d0- 1  for  pure  Pt(lll)  and  Pt(100)  are  given  in  parenthesis. 


Pt„ML/Pt3Co(  111)  (A) 

ptnML/Pt3Co(i00)  (A) 

do-i 

di-2 

d2- 3 

0 

u 

1 

T3 

Pt-Pt 

do-i 

di-2 

d2- 3 

dpt-co 

Pt-Pt 

Pt  1  ML 

2.31 

— 

— 

2.72 

2.75 

1.90 

— 

— 

2.73 

2.81 

Pt2  ML 

2.39 

2.28 

- 

2.72 

2.75 

1.97 

1.94 

- 

2.78 

2.81 

Pt3  ML 

2.39 

2.38 

2.30 

2.72 

2.75 

1.96 

1.96 

1.96 

2.78 

2.81 

Pure  Pt 

2.31 

2.81 

1.93 

2.81 

(2.26)31 

(1.94)48 

atom  underneath  (see  inset  figure).  The  actual  (lOO)-cut  Ll2-Pt3Co 
structure  has  originally  two  Co  atoms  in  an  “alternating”  PtsCo 
substrate’s  layers.  However,  for  consistency  with  the  composition 
of  Co  atoms  with  the  Pt3Co(lll),  then  one  Co  atom  should  be 
moved  to  the  next  lower  substrate  layer  as  shown  by  the  arrows. 
Although  this  leads  to  a  system  with  quite  a  higher  energy,  it  allows 
for  a  suitable  and  reliable  comparison  of  Co  migration  as  a  function 
of  Pt  surface  structure  and  thickness. 

Next,  the  energetics  of  the  migration  of  Co  to  the  surface  for  the 
bare  PtnML/Pt3Co  for  both  (111)  and  (100)  facets  is  obtained  by 
interchanging  the  Pt  (in  higher  layer)  with  the  Co  (in  the  lower 
layer)  and  then  calculating  the  energy  change  (A E)  upon  this  Co 
upshift  (i.e.  migration).  This  scheme  is  shown  in  the  Fig.  3  (inset). 
For  all  systems,  the  reference  energy  for  A E  is  the  total  energy  of  the 
initially  optimized  Pt-skin  structures  as  depicted  in  Figs.  1  and  2. 
The  AE  is  calculated  for  all  values  of  n.  For  the  PtnML/Pt3Co(lll ),  the 
AE  plots  for  Co  migration  to  the  surface  is  shown  in  Fig.  3  and  that 
for  PtnML/Pt3Co(100)  is  given  in  Fig.  4.  First,  for  the  PtnML/Pt3Co(lll ), 
it  can  be  noted  that  there  is  a  very  little  positive  change  (endo¬ 
thermic)  in  migration  energy  (AE  =  0.02-0.06  eV)  when  Co  moves 
within  the  Pt-skin  {indicated  as  Pt  mid-layers  in  the  figure},  while 
there  is  a  large  positive  change  (endothermic)  in  energy 
(AE  =  0.65-0.70  eV)  when  Co  migrates  from  the  subsurface  to  the  Pt- 
skin  surface.  In  general,  this  agrees  with  the  theoretical  study  on 
segregation  of  metals,  which  reports  that  Pt,  and  not  Co,  segregate 
easily  to  the  surface  in  a  Pt-Co  bimetallic  system  [37].  For  the  bare 
PtnML/Pt3Co(100),  the  same  trend  in  Co  migration  to  surface  can  be 
observed.  The  AE  plot  is  shown  in  Fig.  4.  We  note  that  the  AE  for  Co 
migration  in  the  Pt  mid-layers  are  also  minimal,  however,  we  see 
that  there  are  energies  that  are  downhill  (exothermic).  This  is  due 


to  the  fact  that  the  Pt3Co(100)  system  is  formed  from  interchanging 
Co  and  Pt  positions  to  allow  for  the  same  composition  of  Co  and  Pt 
per  layer  as  that  of  Pt3Co(lll )  system,  as  discussed  previously.  Thus, 
when  the  Co  reached  the  so  called  “alternate”  layers  where  they  are 
more  stable,  a  negative  change  in  migration  energy  can  be  ex¬ 
pected.  As  in  PtnML/Pt3Co(lll ),  the  migration  from  the  subsurface  to 
the  Pt  surface  for  all  three  PtnML/Pt3Co(100)  systems  also  entail 
large  energy  cost  (AE  =  0.53-0.72  eV).  Thus,  for  both  (111)  and 
(100)  Pt-skin  structures,  two  generalizations  can  be  deduced:  (1)  a 
minimal  change  in  energy  is  needed  to  move  within  the  Pt-skin  and 
(2)  large  positive  energies  are  required  for  Co  migration  from  sub¬ 
surface  to  the  surface. 

3.2.  Effect  of  oxygen  on  Co  migration 

Because  of  the  suspected  role  of  oxygen  in  the  segregation  of  Co 
to  the  surface  via  formation  of  surface  oxide,  the  adsorption  of 
oxygen  on  all  the  6  Pt-skin  structures  is  conducted.  The  possible 
adsorption  sites  of  0  in  PtnML/Pt3Co(lll)  is  shown  in  Fig.  5.  Bridge 
(b),  hollow  (fee,  hep)  and  top  (t)  sites  are  the  main  symmetric  sites 
and  are  further  split  into  various  other  sites  based  on  the  subsurface 
atoms.  Overall,  five  types  of  bridge  sites  are  identified,  two  for  hep, 
two  for  fee  and  also  two  for  top  sites.  The  adsorption  on  a  (2  x  2) 
surface  unit  cell  leads  to  a  0.25  ML  oxygen  coverage.  The  binding 
energies  are  calculated  based  on  the  following: 

^ads  =  ~  (^slab  +  Eq  atom) 

where  ET  is  the  total  energy  of  the  adsorbed  system,  Eslab  and 
E0  atom  represent  the  total  energies  of  isolated  slab  and  oxygen 


Layer  Number 


Fig.  3.  Co  migration  energy  (AE  in  PtnML/Pt3Co(lll)  systems.  Up  arrows  indicate 
positive  AE  (endothermic)  and  down  arrows  depict  negative  AE  (exothermic).  Sample 
Pt  and  Co  exchange  of  positions  is  illustrated  in  the  inset  figure.  AE  is  obtained  from 
the  total  energy  of  the  new  Co-migrated  PtnML/Pt3Co(lll)  system  referenced  to  the 
initial  PtnML/Pt3Co(lll).  Regions  in  the  plot  where  energies  corresponding  to  Co 
migration  to  the  Pt-skin  surface  and  to  that  within  the  Pt  mid-layers  are  shown  in 
dashed  boxes. 


Fig.  4.  Co  migration  energy  (AE  in  PtnML/Pt3Co(100)  systems.  Up  arrows  indicate 
positive  AE  (endothermic)  and  down  arrows  depict  negative  AE  (exothermic).  AE  is 
obtained  from  the  total  energy  of  the  new  Co-migrated  PtnML/Pt3Co(100)  system 
referenced  to  the  initial  optimized  slab.  The  same  scheme  of  Co  migration  towards  the 
surface  as  that  of  (111)  is  used.  Regions  in  the  plot  where  energies  corresponding  to  Co 
migration  to  the  Pt-skin  surface  and  to  that  within  the  Pt  mid-layers  are  shown  in 
dashed  boxes. 
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(a)  (b)  (c) 


Fig.  5.  (a),  (b)  and  (c)  show  the  top  views  of  PtiML/Pt3Co(lll),  Pt2ML/Pt3Co(lll)  and  Pt3ML/Pt3Co(lll),  respectively.  The  oxygen  atom  adsorption  sites  are  labeled  accordingly.  The 
most  stable  site  is  indicated  in  bold  form. 


atom,  respectively.  Binding  energies  are  given  in  Table  2.  The  most 
stable  oxygen  adsorption  sites  are  indicated  in  bold  form.  It  can  be 
noted  that  the  0  atom  prefers  the  fee  sites,  which  is  the  same  as  the 
most  favorable  adsorption  position  of  oxygen  atom  in  the  pure 
Pt(lll).  However,  the  binding  energies  are  lower  on  the  Pt-skin 
than  on  the  pure  Pt  and  the  Eads  tends  to  increase  upon  increase 
of  the  thickness  of  the  Pt-skin.  The  lower  binding  energy  of  0  atom 
on  other  Pt-3d  transition  metal  alloy  systems  as  compared  to  that  of 
pure  Pt  has  also  been  reported  [38—41  ].  We  further  note  that  the 


Table  2 

Adsorption  sites  and  corresponding  binding  energies  (Eads)  °f  oxygen  atom  on  PtnML/ 
Pt3Co(lll),  PtnML/Pt3Co(100),  Pt(lll)  and  Pt(100).  The  most  stable  adsorption  sites 
are  indicated  in  bold  form.  The  bridge  sites  for  (111)  facet  relaxes  to  the  nearby  fee 
sites  and  the  final  state  is  indicated  in  the  superscript,  b*  for  pure  Pt(lll)  indicates 
the  local  minima  is  found  using  residual  minimization  (RMM-DIIS)  algorithm. 

PtnML/Pt3Co(  111)  PtnML/Pt3Co(  1 00) 


Site  £ads  (eV)  Site  Eads  (eV) 


Pt  1ML 

tl 

-3.76 

tl 

-3.30 

t2 

-2.90 

t2 

-3.32 

hcpl 

-3.72 

hi 

-4.27 

hcp2 

-3.76 

h2 

-4.14 

feel 

-4.12 

bl 

-4.34 

fcc2 

-3.85 

b2 

-4.32 

bl 

-3.85fcc2 

b2 

-3.85fcc2 

b3 

-4.12fccl 

b4 

-4.12fccl 

b5 

-3.85fcc2 

Pt2  ML 

tl 

-2.96 

tl 

-3.44 

t2 

-2.91 

t2 

-3.51 

hcpl 

-3.87 

hi 

-4.31 

hcp2 

-3.89 

h2 

-4.17 

feel 

-4.22 

bl 

-4.44 

fcc2 

-4.30 

b2 

-4.42 

bl 

-4.22fccl 

b2 

-4.22fccl 

b3 

-4.22fccl 

b4 

-4.22fccl 

b5 

-4.22fccl 

Pt3  ML 

tl 

-2.95 

tl 

-3.47 

t2 

-2.98 

t2 

-3.48 

hep 

-3.92 

hi 

-4.15 

feel 

-4.28 

h2 

-4.21 

fcc2 

-4.32 

bl 

-4.48 

bl 

-4.28fccl 

b2 

-4.46 

b2 

-4.28  fccl 

b3 

-4.28  fccl 

b4 

-4.32  fcc2 

b5 

-4.28  fccl 

Pure  Pt 

t 

-3.06 

t 

-3.55 

hep 

-4.01 

h 

-4.23 

fee 

-4.46 

b 

-4.60 

b 

-4.46fcc 

b* 

— 3.74(rmm) 

bridge  sites  are  not  the  global  minima  since  the  oxygen  slides  to 
adjacent  fcc-sites  as  indicated  in  Table  2.  This  is  due  to  the  conju¬ 
gate  gradient  algorithm  used  in  the  optimization  as  mentioned  in 
the  computational  method  section.  This  method  searches  for  the 
global  minimum  and  not  the  local  minima.  This  is  verified  by  the 
test  on  pure  Pt(lll)  using  residual  minimization  method  (RMM- 
DIIS)  [42,43],  known  for  the  local  minima  identification.  The 
binding  energy  for  bridge  site  (b*  in  Table  2)  is  obtained.  The 
location  of  bridge  site  in  Pt(lll)  using  this  algorithm  is  in  agree¬ 
ment  with  previous  work  [44].  Next,  the  slight  destabilization  of  O 
atom  on  catalyst  surface  as  compared  to  the  conventional  Pt  cata¬ 
lyst  is  often  taken  as  a  predictor  of  better  ORR  activity.  The  less 
bound  O  adsorbates  can  facilitate  protonation  leading  to  easier 
formation  of  intermediate  species  in  the  over-all  ORR  reaction 
pathways  [45,46].  Under  0.25  ML  oxygen  coverage,  the  migration 
energy  profile  for  Co  in  PtnML/Pt3Co(lll)  is  now  changed.  In 
Fig.  6(a),  it  is  shown  that  the  Co  migration  in  single  Pt  layer  on 
Pt3Co(lll)  to  the  surface  is  exothermic,  A E  =  -0.03  eV  (while 
endothermic  in  zero-coverage).  The  Co-0  bond  formation  on  the 
surface  is  shown  in  the  inset  figure.  For  the  thicker  Pt-skins  {Pt2- 
3ML/Pt3Co(lll)},  the  AE  are  all  positive.  However,  the  migration 
from  Pt  subsurface  to  the  surface  is  downhill,  that  is,  when 
compared  to  the  previous  step  (see  the  down  arrows  at  the  end  for 
n  =  2  and  n  =  3).  It  is  also  worthwhile  to  note  that  when  the 
migration  of  Co  in  the  zero  oxygen  coverage  (Fig.  3)  and  in  the  0.25 
ML  coverage  (Fig.  6)  are  compared,  the  trends  in  the  movement  of 
Co  within  the  Pt-skin  is  preserved  (same  sequence  of  up  and  down 
arrows).  In  fact,  the  positive  AE  for  the  migration  in  the  middle 


Fig.  6.  Co  migration  energy  (AE  in  PtnML/Pt3Co(lll)  systems  under  0.25  ML  oxygen 
coverage.  Up  arrows  indicate  positive  AE  (endothermic)  and  down  arrows  depict 
negative  AE  (exothermic).  The  same  scheme  of  Co  migration  towards  the  surface  as 
that  of  bare  PtnML/Pt3Co(lll)  slab  is  used.  The  Co-0  bond  formation  on  the  surface  is 
shown  in  the  inset  figures,  corresponding  to  the  Co  migration  from  subsurface  to  the 
surface  (indicated  by  line  arrows). 
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layers  are  greater  in  Pt2-3ML/Pt3(Colll)  under  0.25  ML  oxygen 
coverage  than  in  the  zero  oxygen  coverage.  Thus,  two  important 
characteristics  of  Co  migration  with  varied  Pt-skin  thickness  under 
oxygen  coverage  can  be  derived.  First,  the  thicker  Pt  layers  can 
stabilize  Co  within  the  system,  which  now  confirms  the  experi¬ 
mental  observation  on  stable  PtSkin/Pt3Co  with  2-3  Pt  atomic  layers 
against  Co  dissolution  [12,13].  Secondly,  the  Co  migration  becomes 
downhill  in  energy  only  when  the  subsurface  (one  layer  below  the 
surface)  is  reached.  This  explains  the  experimental  observation  that 
Co  dissolution  is  mostly  prevalent  in  a  1  atom  layer-thick  Pt  skin 
over  Pt-Co  alloy  4,8]. 

Similarly,  the  adsorption  of  O  atom  is  investigated  for  PtnML/ 
Pt3Co(100).  The  adsorption  sites  are  shown  in  Fig.  7.  The  hollow  (h), 
bridge  (b)  and  top  (t)}  sites  are  also  further  split  into  two  differing 
oxygen  positions  (e.g.  bl  and  b2  for  bridge  site,  and  so  on).  The 
binding  energies  are  also  provided  in  Table  2.  The  most  stable  site  is 
labeled  in  bold  form.  It  the  same  way,  the  Eads  in  the  PtnML/ 
Pt3Co(100)  are  lower  than  that  of  the  pure  Pt(100).  The  most  stable 
bridge  site  adsorption  in  Pt(100)  is  also  preserved  in  the  Pt-skin 
systems.  Fig.  8  shows  the  change  in  the  migration  energy,  A E  of 
Co  for  different  Pt  outerlayer  thickness.  For  PtiML/Pt3Co(100),  the 
migration  is  surprisingly  endothermic  with  a  small  positive  energy, 
A E  =  0.04  eV.  In  fact,  the  Co  migration  from  the  subsurface  to  the 
Pt-skin  surface  for  all  thickness  is  endothermic  (shown  by  the  up 
arrows  at  the  end  for  n  =  2  and  n  =  3).  This  is  something  to  be 
looked  into  especially  when  compared  to  the  exothermic  migration 
of  Co  from  Pt  subsurface  to  surface  in  (111)  facet  (down  arrows  at 
the  end  for  n  =  2  and  n  =  3  as  shown  previously  in  Fig.  6).  The 
mechanism  of  this  difference  in  the  Co  migration  with  respect  to 
facet  is  discussed  later.  Next,  the  same  combination  of  exothermic 
and  endothermic  Co  migration  within  the  Pt  mid-layers  for  thicker 
Pt2-3ML/Pt3Co(100)  systems  as  that  with  their  bare  slabs  (Fig.  4)  is 
observed.  Thus,  for  thicker  Pt2-3ML/Pt3Co(100)  systems,  the  oxygen 
also  does  not  affect  the  movement  of  Co  within  the  Pt-skin,  similar 
to  that  of  PtnML/Pt3Co(lll).  In  addition  to  the  presence  of  positive 
AE,  the  above  results  can  explain  the  experimentally  observed 
stability  of  Pt  outerlayer  with  (100)  structure  against  potential 
cycling  [13]. 

3.3.  Mechanism  of  the  effect  of  oxygen 

The  striking  difference  in  the  energetics  of  Co  migration  from 
subsurface  to  Pt-skin  surface  upon  O  adsorption  between  (111)  and 
(100)  surface  structures  necessitates  a  mechanism.  In  this  work, 
charge  density  difference  is  employed  before  and  after  O  adsorp¬ 
tion  for  PtiML/Pt3Co(lll)  and  PtiML/Pt3Co(100)  systems  to  deter¬ 
mine  the  extent  of  electronic  effect  caused  by  oxygen  via  a 


Fig.  8.  Co  migration  energy  (A E  in  PtnML/Pt3Co(100)  systems  under  0.25  ML  oxygen 
coverage.  Up  arrows  indicate  positive  AE  (endothermic)  and  down  arrows  depict 
negative  AE  (exothermic).  The  same  scheme  of  Co  migration  towards  the  surface  as 
that  of  bare  PtnML/Pt3Co(100)  slab  is  used.  The  Co-0  bond  formation  on  the  surface  is 
shown  by  the  inset  figures,  which  correspond  to  Co  migration  from  subsurface  to  the 
surface  (indicated  by  line  arrows). 

visualization  of  charge  flow.  The  charge  difference  is  obtained  using 
the  following  expressions: 

PT  ~  (Po  +  PPt-skin) 

where  pT,  p0,  pPt-Skin  are  the  charge  densities  of  the  O-adsorbed 
PtnML/Pt3Co  systems,  isolated  O  atom  and  the  bare  PtnML/Pt3Co 
systems,  respectively.  Similarly,  for  the  bare  PtnML/Pt3Co  systems, 
the  charge  difference  at  the  first  interlayer  is  obtained  via  the 
following: 

Pt  ~  (Ppt  +  Pa) 

where  pTf  pPt,  pA  are  the  charge  densities  of  the  bare  PtnML/Pt3Co 
systems,  isolated  Pt  layer  and  the  Pt3Co  alloy  slab,  respectively. 
Figs.  9  and  10  show  the  charge  redistribution  plots  for  both  systems. 
First,  it  can  be  noted  that  without  oxygen  coverage  {Figs.  9(a)  and 
10  (a)},  the  charge  accumulates  (dark/red)  within  the  Pt-Co 
bonds,  while,  charge  is  depleted  (light/green)  in  the  Pt  sites.  This 
supports  a  strong  Pt-Co  bonding  in  the  bare  systems,  which  in  turn 
can  account  for  the  large  positive  AE.  Flowever,  when  oxygen  is 
adsorbed  {Figs.  9(b)  and  10(b)},  we  note  that  first,  for  the  O  on  the 
fee  site  ofPtiML/Pt3Co(lll),  the  charge  distribution  is  reversed.  The 
accumulation  of  charge  is  now  observed  in  Pt— O  region,  while 
charges  are  depleted  in  Pt-Co  region  {Fig.  9(b)}.  The  strong  Pt-Co 
bonding  noted  in  the  bare  system  is  lost  under  0.25  ML  oxygen 


Fig.  7.  (a),  (b)  and  (c)  show  the  top  views  of  Pt1ML/Pt3Co(100),  Pt2ML/Pt3Co(100)  and  Pt3ML/Pt3Co(100),  respectively.  The  oxygen  atom  adsorption  sites  are  labeled  accordingly.  The 
most  stable  site  is  indicated  in  bold  form. 
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Depletion  ^ 


Accumulation 


(b)  O  at  fee  site 


Accumulation 


Depletion 


Pt1ML/  Pt3Co(lll)  0/Pt1ML/Pt3Co(lll) 


Fig.  9.  (a)  and  (b)  show  the  charge  redistribution  plot  for  bare  PtiML/Pt3Co(lll)  system  and  O-adsorbed  PtiML/Pt3Co(lll),  respectively.  Darker  regions  (red  online)  indicate 
accumulation  of  charge  and  lighter  regions  (green  online)  show  depletion  of  charges.  The  site  of  the  oxygen  adsorption  is  shown.  (For  interpretation  of  the  references  to  color  in  this 
figure  legend,  the  reader  is  referred  to  the  web  version  of  this  article.) 


coverage.  In  the  case  of  PtiML/Pt3Co(100),  as  depicted  in  Fig.  10, 
although  the  same  charge  redistribution  as  that  of  PtiML/Pt3Co(lll ) 
can  be  observed  (i.e.  accumulation  in  Pt-0  and  depletion  in  Pt-Co 
region),  the  bridge  site  adsorption  of  oxygen  allows  for  a  depletion 
of  charge  only  along  two  Pt-Co  bonds,  leaving  the  charge  distri¬ 
bution  in  the  other  Pt-Co  bonds  unchanged.  The  fee  site  adsorption 
of  0  atom  in  the  PtiML/Pt3Co(lll),  on  the  other  hand,  affects  all  the 
Pt-Co  bonds.  This  contribute  to  the  positive  A E  for  Co  migration  in 
(100)  facet,  and  a  negative  A E  for  (111)  facet.  It  has  to  be  noted, 
however,  that  while  at  (100),  the  Co  migration  is  endothermic,  the 
magnitude  of  the  AE  is  minimal.  Finally,  we  further  note  that  by 
virtue  of  the  charge  difference  expressions  above,  the  Pt-Co  bonds 
deeper  into  the  bulk  where  the  charge  densities  rarely  change,  do 
not  show  charge  density  accumulation  (or  depletion)  due  to 
canceling  out  of  same  charge  densities. 

3.4.  Effect  of  Pt-skin  composition 

In  the  quest  for  further  lowering  Pt  loading,  a  possible  substitute 
metal  for  the  thicker  Pt-skin/Pt3Co  systems  is  investigated.  In 
Fig.  11(a),  we  show  the  new  Pt2ML/Pt3Co(lll)  system  with  MBL- 


substituted  Pt  mid-layer.  The  system  is  denoted  now  as  PtiML / 
MBL/Pt3Co(lll ),  where  MBL  =  Ru,  Rh,  Pd,  Os,  Ir.  The  (111 )  surface  is 
used  due  to  the  stable  hexagonal  structure  of  MBLs.  The  insertion  of 
the  MBL  leads  to  a  more  stable  PtiML/MBL/Pt3Co(lll)  than  the 
segregated  MBL-counterpart:  MBL/PtiML/Pt3Co(lll).  This  is  evi¬ 
denced  by  the  lower  total  energy  of  the  former  system  (EA)  as 
shown  in  Table  3  and  furthermore  by  the  relative  stability  as  ob¬ 
tained  by  the  following  expression: 

A  Es  =  Ea-Eb 

Flere,  the  EB  is  the  total  energy  of  the  segregated  counterpart, 
MBL/PtiMi/Pt3Co(lll).  The  same  scheme  for  Co  migration  to  sur¬ 
face  as  in  Pt2ML/Pt3Co(lll)  is  used.  The  AE  plots  under  0.25  ML 
oxygen  coverage  are  shown  in  Fig.  11(b).  The  five  PtiML/MBL/ 
Pt3Co(lll)  systems  are  compared  with  that  of  Pt2ML/Pt3Co(lll).  In 
general,  the  migration  of  Co  in  the  Pt  mid-layer  is  endothermic  for 
all  systems.  Next,  the  migration  from  the  Pt  mid-layer  to  the  surface 
is  exothermic  (referenced  to  the  previous  step),  showing  the  same 
Co  migration  behavior  as  in  the  Pt2ML/Pt3Co(lll).  With  respect  to 
zero  energy  reference,  however,  certain  MBL-substituted  Pt-skin 
systems  exhibit  higher  positive  migration  energies  than  the  Pt2Mi7 


Depletion 

Accumulation 


Accumulation 


O  at  bridge 


Accumulation 


Depletion 


PW  Pt3Co(100)  O/PW  Pt3Co(100) 


Fig.  10.  (a)  and  (b)  show  the  charge  redistribution  plot  for  bare  PtiML/Pt3Co(100)  system  and  O-adsorbed  PtiML/Pt3Co(100),  respectively.  Darker  regions  (red  online)  indicate 
accumulation  of  charge  and  lighter  regions  (green  online)  show  depletion  of  charges.  The  site  of  the  oxygen  adsorption  is  shown.  (For  interpretation  of  the  references  to  color  in  this 
figure  legend,  the  reader  is  referred  to  the  web  version  of  this  article.) 
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(a) 


Pt1Ml/MBL/Pt3Co(lU) 
MBL=Rh,  Ru,  Os,  Pd,  Ir 


(b) 


Fig.  11.  (a)  Migration  Barrier  Layer  (MBL)-substituted  Pt-skin  system:  PtiML/MBL/Pt3Co(lll),  where  MBL  =  Ru,  Rh,  Pd,  Os,  Ir  and  (b)  Co  migration  energy  (A E  plots  for  PtiMiVMBL/ 
Pt3Co(lll)  and  Pt2ML/Pt3Co(lll)  under  0.25  ML  oxygen  coverage. 


Pt3Co(lll).  These  are  when  the  MBLs  are  made  of  Os  and  Pd  while 
those  that  exhibits  negative  change  in  migration  energies  are 
composed  of  MBL  =  Ru,  Rh,  Ir.  Among  these,  Ir  renders  the  least 
negative  change  in  energy.  Overall,  based  on  the  stability  against  Co 
migration  to  surface,  the  PtiML/Os/Pt3Co(lll)  poses  as  the  most 
stable  system.  Furthermore,  as  pointed  out  previously,  the  binding 
energy  of  O  atom  can  serve  as  good  predictor  of  ORR  activity  in  the 
Pt-skin  systems  [45,46].  Hence,  the  binding  energies  of  the  0  atom 
in  these  systems  are  also  calculated  and  are  shown  in  Table  3.  It  can 
be  very  well  noted  that  the  O  adsorption  are  less  stable  in  all  the 
MBL-substituted  Pt-skin  as  compared  to  that  of  Pt2ML/Pt3Co(lll) 
system.  Therefore,  PtiML/MBL/Pt3Co(lll)  where  MBL  =  Os,  Pd 
render  a  more  stable  and  potentially  more  reactive  Pt-skin  ORR 
catalysts. 

3.5.  Stability  against  Pt  dissolution 

Apart  from  the  Co  migration  to  the  surface  with  Co-0  bond 
formation,  the  dissolution  of  Pt  in  both  the  thicker  Pt-skins  and  the 
MBL-substituted  Pt-skin  systems  is  tested.  In  the  kinetic  model  of  Pt 
dissolution  in  PEMFC  [47,48],  three  possible  reactions  have  been 
well-accepted: 

Pt  Pt2+  +  2e~  (Pt  dissolution)  ( 1 ) 

Pt  +  H2O^PtO  +  2H+  +  2e_  (Pt  oxide  film  formation)  (2) 

PtO  +  2H+  -►  Pt2+  +  H20  (chemical  dissolution  of  Pt  oxide) 

(3) 


Table  3 

Stability  data  for  PtiML/MBL/Pt3Co(lll)  systems  with  respect  to  MBL/PtiML/ 
Pt3Co(lll)  systems.  EA  is  the  total  energy  of  the  former  and  EB  is  the  total  energy  of 
the  latter.  AE  (eV)  is  EA  -  EB. 


A 

Ea 

B 

Eb 

A E  (eV) 

PtiML/Ir/Pt3Co(lll) 

-154.76 

Ir/Pt1ML/Pt3Co(  111) 

-153.04 

-1.72 

Pt1ML/Os/Pt3Co(lll) 

-163.66 

Os/Pq  ML/Pt3Co(  ill) 

-161.03 

-2.63 

Pt1ML/Pd/Pt3Co(lll) 

-141.12 

Pd/Pt1ML/Pt3Co(l  1 1 ) 

-140.91 

-0.20 

PtiML/Rh/Pt3Co(lll) 

-149.24 

Rh/Pq  ML/Pt3Co(  111) 

-147.63 

-1.61 

Pt1ML/Ru/Pt3Co(lll) 

-156.39 

Ru/Pq  ML/Pt3Co(  111) 

-153.82 

-2.57 

Therefore,  there  are  two  ways  by  which  Pt  is  lost  to  the  solution, 
the  first  is  the  electrochemical  dissolution  of  Pt  according  to  reac¬ 
tion  (1)  and  the  second  is  the  chemical  dissolution  of  PtO  film  ac¬ 
cording  to  reaction  (3).  The  fact  that  (3)  is  regarded  as  a  slow 
reaction  and  that  it  occurs  at  high  potentials  (when  the  Pt  surface  is 
covered  with  an  oxide  film),  the  viability  of  reaction  (1 )  as  basis  for 
evaluation  of  Pt  dissolution  at  zero  and  low  oxygen  coverage  (0.25 
ML)  is  clear.  When  the  dissolution  of  Pt  is  referenced  to  the  stan¬ 
dard  hydrogen  electrode  (SHE),  then  the  redox  reaction  (1)  {with 
full  description  of  net  reaction  for  the  reference  and  working 
electrode  system}  can  be  rewritten  as: 

Ptn+  +  n/2H2  -►  Pt  +  nH+  (4) 

and  the  Gibbs  free  energy  change  of  the  reaction  is: 

AG  =  /iPt  +  Mn H+  -  Pptn+  -  Mn/ 2H2  (3) 

where  ^Pt,  /iPtn+  are  the  chemical  potentials  of  pure  metal  Pt  and  Pt 
ions  (and  that  equilibrium  exists  between  Pt  in  the  system  and  the 
dissolved  Pt  in  the  solution)  and  ^nH+,  gn/2u2  are  the  chemical  po¬ 
tentials  of  H+  and  H2  (at  standard  conditions).  The  AG  is  related  to 
the  electrochemical  potential  by: 


AG  =  neUpt  =  Mpt  +  M„h+  -  MPt"+  -  Mn/2H2  (6) 

For  other  type  of  Pt  systems,  denoted  as  Pt’,  equation  (6)  is  thus: 

neUpt  =  MPt-  +  MnH+  -  Mptn+  -  Mn/2H2  (7) 

and  the  shift  in  the  electrode  potential  is  (after  cancellation  of  same 
terms): 


ne(Lfpt’  U pt)  —  MPt  Mpt  (3) 

In  this  work,  the  chemical  potentials  are  calculated  based  on  Pt 
vacancy  formation  energy,  £des,  in  the  presence  of  0.25  ML  oxygen 
coverage: 


^des  — 


^Ept  +  E( 


0/Pt_vac 


-Et 


O/Pt 


(9) 
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(a)  (b)  (c) 


Pt2ML/ Pt3Co(lll)  Pt3ML/Pt3Co(lll)  Pt1Ml/MBL/Pt3Co(lll) 

MBL  =  Os,  Ir,  Pd 


Fig.  12.  (a),  (b)  and  (c)  show  the  atomic  structure  of  Pt2ML/Pt3Co(lll),  Pt3ML/Pt3Co(lll)  and  PtiML/M/Pt3Co(lll)  with  1  Pt  atom  vacancy,  respectively.  The  Pt,  Co  and  0  atoms  are 
shown. 


where  EPt,  Eo/Pt_vac  and  E0/Pt  are  the  total  energies  of  isolated  Pt 
atom,  O/Pt  with  1  Pt  vacancy  and  the  O/Pt  slab  systems,  respec¬ 
tively.  Fig.  12(a)-(c)  shows  the  optimized  Pt2-3ML/Pt3Co(lll)  and 
PtiML/MBL/Pt3Co(lll )  systems  with  1  Pt  atom  vacancy.  Similarly  for 
other  types  of  Pt  systems,  Pt : 

^des  =  (Ept  +  £o/Pt’_vac)  -  £0/Pt 
Thus,  finally: 

Mpt’  -  Mpt  =  ^des  _  ^des  =  (^O/Pt’.vac  -  ^O/Pt  )  _  (^O/Pt.vac  ~  ^O/Pt) 

(ID 

MJ  =  Upt  -  UPt  =  (Mpt  -  MPt )/ne  (12) 

Table  4  shows  the  chemical  potential  and  electrode  potential 
shifts  for  several  Pt  systems.  The  reference  system  is  Pt(lll).  The 
thicker  Pt-skin  systems  {Pt2-3ML/Pt3Co(lll)}  have  a  very  minimal 
negative  shift  in  the  electrode  potential  (i.e.  -0.09  V  and  -0.06  V). 
Thus,  it  can  be  considered  that  their  stability  against  Pt  dissolution 
should  almost  be  the  same  as  that  of  pure  Pt(lll).  For  the  PtiML / 
MBL/Pt3Co(lll)  systems,  where  MBL  =  Os,  Pd  and  Ir,  Pt  dissolution 
occurs  at  higher  potentials  as  compared  to  that  of  Pt2-3Mi7 
Pt3Co(lll).  Interestingly,  those  composed  of  Os  and  Pd  are  in  fact 
exhibiting  a  positive  potential  shift.  This  suggests  better  stability 
than  the  pure  Pt(lll)  itself.  The  PtiML/Os/Pt3Co(lll)  is  the  most 
stable  against  Pt  dissolution  with  a  high  positive  potential  shift  of 
0.16  V.  Since  Os  is  four  to  five  times  cheaper  than  Pt,  then  this  MBL 
offers  a  potential  low-cost  Pt-skin  composition  for  a  Pt-skin  -based 
ORR  catalyst. 

4.  Summary 

First-principles  study  based  on  DFT  is  conducted  to  evaluate  the 
stability  of  Pt-  skin/Pt3Co  systems  against  Co  migration  and  Pt 
dissolution  as  a  function  of  Pt-skin  thickness,  surface  structure  and 


Table  4 

Adsorption  energies  of  0  atom  on  the  Pt  systems;  and  the  chemical  potential  and 
electrochemical  potential  shift  for  Pt  systems. 


Pt’  Systems 

£ads 

PPt  ~  Pet 

A  U 

Pt2ML/Pt3Co(  111) 

-4.30 

-0.18 

-0.09 

Pt3ML/Pt3Co(  111) 

-4.32 

-0.11 

-0.06 

Pt1ML/Os/Pt3Co(lll) 

-4.20 

0.31 

0.16 

Pt1ML/Pd/Pt3Co(  111) 

-4.27 

0.06 

0.03 

PtiML/Ir/Pt3Co(  111) 

-4.26 

-0.08 

-0.04 

composition.  It  is  found  that  the  thicker  Pt-skin  (Pt2-3ML/Pt3Co)  in 
both  the  (100)  and  (111)  facets  under  0.25  ML  oxygen  coverage  are 
more  stable  than  PtiML/Pt3Co  counterpart,  confirming  the  observed 
stability  in  thicker  Pt-skin  systems  on  Pt-Co  alloys  in  experiments. 
The  enhanced  stability  is  attributed  to  the  positive  change  in  energy 
(A E  for  Co  migration  within  the  Pt  mid-layers  as  well  as  the  no 
found  effect  of  adsorbed  O  atom  on  the  Co  movement  within  the  Pt- 
skin.  Flowever,  it  has  to  be  noted  that  the  behavior  of  migration  of 
Co  from  the  subsurface  to  the  Pt-skin  surface  is  different  for 
different  facets.  We  note  that  the  Co  migration  is  exothermic  in 
(111)  and  endothermic  in  (100)  surface  structures.  Charge  differ¬ 
ence  plots  reveal  that  upon  oxygen  adsorption,  the  charge  deple¬ 
tion  in  Pt-Co  bond  is  less  in  (100)  facet  compared  to  that  of  (111). 
This  is  due  to  the  bridge  site  adsorption  of  oxygen  in  PtiML/ 
Pt3Co(100)  as  compared  to  the  hollow  site  (fee)  adsorption  in  PtiML/ 
Pt3Co(lll).  Furthermore,  the  stability  against  Pt  dissolution  in  the 
thicker  Pt-skin  system,  Pt2-3ML/Pt3Co(lll),  is  evaluated.  It  is  noted 
that  a  very  minimal  negative  electrochemical  potential  shift  is 
exhibited  by  the  Pt-skin  systems,  signifying  a  stability  close  to  the 
reference  system,  Pt(lll).  Lastly,  the  Pt-skin  composition  is  varied 
by  substituting  Migration  Barrier  Layer  (MBL)  for  the  Pt  mid-layer 
in  the  thicker  Pt-skin  {i.e.  PtiML/MBL/Pt3Co(lll)}.  The  MBL  is  a  4d 
or  a  5d  transition  metal:  Ru,  Rh,  Ir,  Os,  Pd.  The  stability  against  Co 
migration  and  Pt  dissolution  under  oxygen  coverage  is  also  tested. 
It  is  observed  that  on  all  the  PtiML/MBL/Pt3Co(lll)  systems,  the 
oxygen  binding  energy  is  less  than  on  the  corresponding  Pt2ML/ 
Pt3Co(lll)  as  well  as  on  the  Pt(lll).  The  Co  migration  to  surface  in 
PtiML/MBL/Pt3Co(lll),  where  M  is  Os  and  Pd  is  best  inhibited  as 
compared  to  that  of  the  Pt2ML/Pt3Co(lll)  and  the  other  MBL- 
substituted  Pt-skins,  due  to  their  higher  positive  magnitude  for 
A E.  The  electrochemical  potential  shift  for  Pt  dissolution  for  the 
thicker  Pt-skins  {Pt2-3ML/Pt3Co(lll)}  and  for  MBL-substituted  Pt 
skins  {PtiML/MBL/Pt3Co(lll)}  shows  that  the  thicker  Pt-skins  have 
almost  the  same  Pt  stability  as  that  of  Pt( 111 ).  It  is  noteworthy  that 
the  MBL-substituted  Pt-skin  systems  are  in  fact  more  stable  than 
the  thicker  Pt-skin  systems  against  Pt  dissolution.  Screening  of 
MBLs  based  on  the  electrochemical  potential  shift  for  Pt  dissolution 
picks  out  M  =  Os  and  Pd  as  the  best  substitution  metal  due  to  the 
positive  electrode  potential  shift.  This  means  that  Pti  ml/Os/ 
Pt3Co(lll)  and  PtiML/Pd/Pt3Co(lll)  are  more  stable  than  the 
Pt(lll)  itself.  Strikingly,  PtiML/Os/Pt3Co(lll)  gives  highest  positive 
potential  shift  of  ~0.16  V.  Because  Os  is  much  cheaper  than  Pt  and 
that  it  gives  both  the  highest  positive  electrode  potential  shift  for  Pt 
dissolution  and  the  largest  magnitude  for  AE  for  Co  migration  than 
the  rest  of  the  MBLs,  then  Os  poses  a  highly  viable  compositional 
material  for  a  low-cost,  stable  and  active  Pt-skin  ORR  catalyst  for 
PEMFC. 
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